Structural basis for the mechanism of interaction of SARS-CoV-2 B.1.640.2 variant RBD with the host receptors hACE2 and GRP78.

The inflicted chaos instigated by the SARS-CoV-2 (Severe acute respiratory syndrome coronavirus 2) globally continues with the emergence of novel variants. The current global outbreak is aggravated by the manifestation of novel variants, which affect the effectiveness of the vaccine, attachment with hACE2 (human Angiotensin-converting enzyme 2) and immune evasion. Recently, a new variant named University Hospital Institute (IHU) (B.1.640.2) was reported in France in November 2021 and is spreading globally affecting public healthcare. The B.1.640.2 SARS-CoV-2 strain revealed 14 mutations and 9 deletions in spike protein. Thus, it is important to understand how these variations in the spike protein impact the communication with the host. A protein coupling approach along with molecular simulation protocols was used to interpret the variation in the binding of the wild type (WT) and B.1.640.2 variant with hACE2 and Glucose-regulating protein 78 (GRP78) receptors. The initial docking scores revealed a stronger binding of the B.1.640.2-RBD with both the hACE2 and GRP78. To further understand the crucial dynamic changes, we looked at the structural and dynamic characteristics and also explored the variations in the bonding networks between the WT and B.1.640.2-RBD (receptor-binding domain) in association with hACE2 and GRP78, respectively. Our findings revealed that the variant complex demonstrated distinct dynamic properties in contrast to the wild type due to the acquired mutations. Finally, to provide conclusive evidence on the higher binding by the B.1.640.2 variant the TBE was computed for each complex. For the WT with hACE2 the TBE was quantified to be-61.38 ± 0.96 kcal/mol and for B.1.640.2 variant the TBE was estimated to be -70.47 ± 1.00 kcal/mol. For the WT-RBD-GRP78 the TBE -was computed to be 32.32 ± 0.56 kcal/mol and for the B.1.640.2-RBD a TBE of -50.39 ± 0.88 kcal/mol was reported. This show that these mutations are the basis for higher binding and infectivity produced by B.1.640.2 variant and can be targeted for drug designing against it.Communicated by Ramaswamy H. Sarma.

Plant-derived active compounds as a potential nucleocapsid protein inhibitor of SARS-CoV-2: an in-silico study.

The coronavirus disease 2019 (COVID-19) is caused by SARS-CoV-2. This virus has a high mismatch repair proofreading ability due to its unique exonuclease activity, making it knotty to treat. The nucleocapsid protein can serve as a potential antiviral drug target, as this protein is responsible for multiple captious functions during the viral life cycle. Herein, we have investigated the potential to repurpose active antiviral compounds of plant origins for treating the SARS-CoV-2 infection. In the present study, we followed the molecular docking methodology to screen druggable natural plants' active compounds against the nucleocapsid protein of SARS-CoV-2. The virtual screening of all 68 compounds revealed that the top seven active compounds, such as withanolide D, hypericin, silymarin, oxyacanthine, withaferin A, Acetyl aleuritolic acid, and rhein, exhibit good binding affinity with druggable ADME properties, toxicity, and Pass prediction. The stability of the docked complexes was studied by conducting molecular simulations of 100 ns. MM-GBSA calculated the binding free energy uncovered that withanolide D, hypericin, and silymarin result in highly stable binding conformations in three different sites of the nucleocapsid protein. However, further investigation is needed in order to validate the candidacy of these inhibitors for clinical trials. HighlightsNatural plants' active compounds may aid in the inhibition of SARS-CoV-2 replication and COVID-19 therapeutics.Hypericin, silymarin, withanolide D, oxyacanthine, withaferin A, Acetyl aleuritolic acid, and rhein are effective against SARS-CoV-2 N protein.Studied natural plants' active compounds could be useful against COVID-19 and its associated organs comorbidities.ADMET properties of selected compounds favor these compounds as druggable candidates.Communicated by Ramaswamy H. Sarma.

Bioinformatic Screening of Compounds from Iranian Lamiaceae Family Members against SARS-CoV-2 Spike Protein

Background: COVID-19 (coronavirus disease 2019) is still a major challenge worldwide. The disease is caused by binding the coronavirus to ACE2 receptors on lung cells, infecting the cells and triggering the onset of symptoms. The prevention of such a binding in which the virus is eventually unable to enter the cell could be a promising therapeutic approach.Methods: In this in silico study, 306 compounds of Lamiaceae family native in Iran (native Mints) were retrieved from several databases as 3D structures, and after that molecular docking and virtual screening, the compounds with inhibitory potential were selected in terms of free energy binding against the spike protein of the virus. The pharmacokinetic profile of selected compounds was evaluated, and by molecular dynamic simulation and MM/PBSA, four compounds were further assessed for binding affinities against the receptor-binding domain of the spike.Results: The results showed the Catechin gallate and Perovskone B from Stachys and Salvia genus generated a stronger binding affinity, and therefore could act as potential inhibitory compounds of RBD of the SARS-CoV-2 spike protein.Conclusion: This study revealed that some members of the Lamiaceae family could be employed to inhibit SARS-CoV-2 activity through interaction with spike protein and therefore could be used for further investigation in vitro and in vivo.

Towards a free energy-based elastic network model and its application to the SARS-COV2 binding to ACE2.

Classical normal mode analysis (cNMA) is a standard method for studying the equilibrium vibrations of macromolecules. A major limitation of cNMA is that it requires a cumbersome step of energy minimization that also alters the input structure significantly. Variants of normal mode analysis (NMA) exist that perform NMA directly on PDB structures without energy minimization, while maintaining most of the accuracy of cNMA. Spring-based NMA (sbNMA) is such a model. sbNMA uses an all-atom force field as cNMA does, which includes bonded terms such as bond stretching, bond angle bending, torsional, improper, and non-bonded terms such as van der Waals interactions. Electrostatics was not included in sbNMA because it introduced negative spring constants. In this work, we present a way to incorporate most of the electrostatic contributions in normal mode computations, which marks another significant step toward a free-energy-based elastic network model (ENM) for NMA. The vast majority of ENMs are entropy models. One significance of having a free energy-based model for NMA is that it allows one to study the contributions of both entropy and enthalpy. As an application, we apply this model to study the binding stability between SARS-COV2 and angiotensin converting enzyme 2 (or ACE2). Our results show that the stability at the binding interface is contributed nearly equally by hydrophobic interactions and hydrogen bonds.

PACAP key interactions with PAC1, VPAC1, and VPAC2 identified by molecular dynamics simulations.

The neuropeptide pituitary adenylate cyclase-activating polypeptide (PACAP) belongs to the glucagon/secretin family. PACAP interacts with the pituitary adenylate cyclase-activating polypeptide receptor type 1 (PAC1) and vasoactive intestinal peptide receptors 1 and 2 (VPAC1 and VPAC2), exhibiting functions in the immune, endocrine, and nervous systems. This peptide is upregulated in numerous instances of brain injury, acting as a neuroprotective agent. It can also suppress HIV-1 and SARS-CoV-2 viral replication in vitro. This work aimed to identify, in each peptide-receptor system, the most relevant residues for complex stability and interaction energy communication via Molecular Dynamics (MD), Free Energy calculations, and Protein-energy networks, thus revealing in detail the underlying mechanisms of activation of these receptors. Hydrogen bond formation, interaction energies, and computational alanine scanning between PACAP and its receptors showed that His1, Asp3, Arg12, Arg14, and Lys15 are crucial to the peptide's stability. Furthermore, several PACAP interactions with structurally conserved positions deemed necessary in GPCR B1 activation, including Arg2.60, Lys2.67, and Glu7.42, were significant for the peptide's stability within the receptors. According to the protein-energy network, the connection between Asp3 of PACAP and the receptors' conserved Arg2.60 represents a critical energy communication hub in all complexes. Additionally, the ECDs of the receptors were also found to function as energy communication hubs for PACAP. Although the overall binding mode of PACAP in the three receptors was found to be highly conserved, Arg12 and Tyr13 of PACAP were more prominent in complex with PAC1, while Ser2 of PACAP was with VPAC2. The detailed analyses performed in this work pave the way for using PACAP and its receptors as therapeutic targets.Communicated by Ramaswamy H. Sarma.

AI-Driven De Novo Design and Molecular Modeling for Discovery of Small-Molecule Compounds as Potential Drug Candidates Targeting SARS-CoV-2 Main Protease.

Over the past three years, significant progress has been made in the development of novel promising drug candidates against COVID-19. However, SARS-CoV-2 mutations resulting in the emergence of new viral strains that can be resistant to the drugs used currently in the clinic necessitate the development of novel potent and broad therapeutic agents targeting different vulnerable spots of the viral proteins. In this study, two deep learning generative models were developed and used in combination with molecular modeling tools for de novo design of small molecule compounds that can inhibit the catalytic activity of SARS-CoV-2 main protease (Mpro), an enzyme critically important for mediating viral replication and transcription. As a result, the seven best scoring compounds that exhibited low values of binding free energy comparable with those calculated for two potent inhibitors of Mpro, via the same computational protocol, were selected as the most probable inhibitors of the enzyme catalytic site. In light of the data obtained, the identified compounds are assumed to present promising scaffolds for the development of new potent and broad-spectrum drugs inhibiting SARS-CoV-2 Mpro, an attractive therapeutic target for anti-COVID-19 agents.

Marine drugs as putative inhibitors against non-structural proteins of SARS-CoV-2: an in silico study.

INTRODUCTION: Coronavirus disease 2019 (COVID-19) is an unprecedented pandemic, threatening human health worldwide. The need to produce novel small-molecule inhibitors against the ongoing pandemic has resulted in the use of drugs such as chloroquine, azithromycin, dexamethasone, favipiravir, ribavirin, remdesivir and azithromycin. Moreover, the reports of the clinical trials of these drugs proved to produce detrimental effects on patients with side effects like nephrotoxicity, retinopathy, cardiotoxicity and cardiomyopathy. Recognizing the need for effective and non-harmful therapeutic candidates to combat COVID-19, we aimed to develop promising drugs against SARS-COV-2. DISCUSSION: In the current investigation, high-throughput virtual screening was performed using the Comprehensive Marine Natural Products Database against five non-structural proteins: Nsp3, Nsp5, Nsp12, Nsp13 and Nsp15. Furthermore, standard precision (SP) docking, extra precision (XP) docking, binding free energy calculation and absorption, distribution, metabolism, excretion and toxicity studies were performed using the SchrÓ§dinger suite. The top-ranked 5 hits obtained by computational studies exhibited to possess a greater binding affinity with the selected non-structural proteins. Amongst the five hits, CMNPD5804, CMNPD20924 and CMNPD1598 hits were utilized to design a novel molecule (D) that has the capability of interacting with all the key residues in the pocket of the selected non-structural proteins. Furthermore, 200 ns of molecular dynamics simulation studies provided insight into the binding modes of D within the catalytic pocket of selected proteins. CONCLUSION: Hence, it is concluded that compound D could be a promising inhibitor against these non-structural proteins. Nevertheless, there is still a need to conduct in vitro and in vivo studies to support our findings.

Molecular insights to the binding interactions of APNS containing HIV-protease inhibitors against SARS-CoV-2 Mpro: an in silico approach towards drug repurposing.

SARS-CoV-2 Mpro is one of the most vital enzymes of the new coronavirus-2 (SARS-CoV-2) and is a crucial target for drug discovery. Unfortunately, there is not any potential drugs available to combat the action of SARS-CoV-2 Mpro. Based on the reports HIV-protease inhibitors can be applied against the SARS by targeting the SARS-CoV-1 Mpro, we have chosen few clinically trialed experimental and allophenylnorstatine (APNS) containing HIV-protease inhibitors (JE-2147, JE-533, KNI-227, KNI-272 & KNI-1931), to examine their binding affinities with SARS-CoV-2 Mpro and to assess their potential to check for a possible drug candidate against the protease. Here, we have chosen a methodology to understand the binding mechanism of these five inhibitors to SARS-CoV-2 Mpro by merging molecular docking, molecular dynamics (MD) simulation and MM-PBSA based free energy calculations. Our estimations disclose that JE-2147 is highly effective (ΔGBind = -28.31 kcal/mol) due to an increased favorable van der Waals (ΔEvdw) interactions and decreased solvation (ΔGsolv) energies between the inhibitor and viral protease. JE-2147 shows a higher level of interactions as compared to JE-533 (-6.85 kcal/mol), KNI-227 (-18.36 kcal/mol), KNI-272 (-15.69 kcal/mol) and KNI-1931 (-21.59 kcal/mol) against SARS-CoV-2 Mpro. Binding contributions of important residues (His41, Met49, Cys145, His164, Met165, Glu166, Pro168, Gln189, etc.) from the active site or near the active site regions with ≥1.0 kcal/mol suggest a potent binding of the inhibitors. It is anticipated that the current study of binding interactions of these APNS containing inhibitors can pitch some valuable insights to design the significantly effective anti-SARS-CoV-2 Mpro drugs.Communicated by Ramaswamy H. Sarma.

Computational Investigations for Identification of Bioactive Molecules from Baccaurea ramiflora and Bergenia ciliata as Inhibitors of SARS-CoV-2 Mpro

In this study, a hybrid compound library of 72 phytocompounds from two antiviral medicinal plants (Baccaurea ramiflora and Bergenia ciliata) was computationally investigated for their inhibitory potential against SARS-CoV-2 Mpro. Molecular docking showed that 6-O-vanilloylicariside B5, 6-O-vanilloylisotachioside, leucoanthocyanidin 4-(2-galloyl), and p-hydroxybenzoyl bergenin has good binding affinity for Mpro. However, p-hydroxybenzoyl bergenin did not bind at the catalytic cavity. The RMSD and RMSF data obtained from 100 ns MD simulations revealed stable protein–ligand complexes for 6-O-vanilloylicariside B5, 6-O-vanilloylisotachioside, leucoanthocyanidin 4-(2-galloyl). Ligand trajectory study found 6-O-vanilloylisotachioside and leucoanthocyanidin 4-(2-galloyl) to be stable. Studies on ligand interaction profile and timeline interaction profile showed that 6-O-vanilloylisotachioside and leucoanthocyanidin 4-(2-galloyl) interacted with HIS41–CYS145 dyad and other crucial amino acids of the catalytic site cavity during the entire 100 ns MD simulations. Molecular mechanics generalized born solvent accessibility binding free energy calculations, density functional theory analysis, quantitative structure–property relationship studies, and ADMET profiling of 6-O-vanilloylisotachioside and leucoanthocyanidin 4-(2-galloyl) supported the results generated by molecular docking and MD simulations studies. Based on the current computational investigations, we conclude that that 6-O-vanilloylisotachioside of B. ramiflora and leucoanthocyanidin 4-(2-galloyl) of B. ciliata are two potential inhibitors of SARS-CoV-2 Mpro that are worthy of further investigations.

Investigate the designing of eutectic mixture using DFT computation and evaluate their biological potency using molecular docking

Researchers are exploring the eutectic mixture because of their obvious great potential in various disciplines. Herein, authors have presented the DFT calculations, molecular docking and QSAR results for designed eutectic mixtures (EMs) using thiourea and resorcinol on taking different equivalent ratio. Authors have used Jakob et al. method to determine the melting temperature of the systems or EMs theoretically. Thermodynamic parameteres such as the free energy, enthalpy, and other energy of the EMs at room temperature are determined through DFT calculations using Gaussian. Authors have also calculated the physiochemical descriptors of various eutectic mixture based on DFT calculations. Further, molecular docking of the designed EMs is carried out to investigate their biological potential for inhibition of the Mpro of SARS-CoV-2. © 2023 Elsevier B.V.

An Evaluation of ASEAN Renewable Energy Path to Carbon Neutrality

Maintaining our standard of living and keeping the economy running smoothly is heavily reliant on a consistent supply of energy. Renewable energy systems create abundant energy by utilizing resources such as the sun, wind, earth, and plants. The demand for renewable energy is increasing, despite power scarcity, pollution, and climate change posing challenges to long-term development in the Association of Southeast Asian Nations (ASEAN), which has seen significant social and economic growth in recent years. To achieve its 23% renewable energy (RE) target, ASEAN can develop solar photovoltaic (PV) electricity. Members of the ASEAN have established regulations and incentives to encourage individuals and businesses to use renewable energy in the future. This paper explores Southeast Asian countries' comprehensive fossil-free energy options, the region's renewable energy potential, current capacity, goals, and energy needs. Through the ASEAN Plan of Action for Energy Cooperation (APAEC) 2016–2025 and the ASEAN Declaration on Renewable Energy, ASEAN is committed to reducing its greenhouse gas emissions and promoting sustainable development aligning with the Paris Agreement's aim to limit global warming to well below 2 degrees Celsius above pre-industrial levels. Results show that decarbonizing the region's energy system is possible, but current policies and actions must be altered to reach that target level. Further research is necessary to optimize the ASEAN region's renewable resource technical potential and commercial viability with available technology.

Template Entrance Channel as Possible Allosteric Inhibition and Resistance Site for Quinolines Tricyclic Derivatives in RNA Dependent RNA Polymerase of Bovine Viral Diarrhea Virus.

The development of potent non-nucleoside inhibitors (NNIs) could be an alternate strategy to combating infectious bovine viral diarrhea virus (BVDV), other than the traditional vaccination. RNA-dependent RNA polymerase (RdRp) is an essential enzyme for viral replication; therefore, it is one of the primary targets for countermeasures against infectious diseases. The reported NNIs, belonging to the classes of quinolines (2h: imidazo[4,5-g]quinolines and 5m: pyrido[2,3-g] quinoxalines), displayed activity in cell-based and enzyme-based assays. Nevertheless, the RdRp binding site and microscopic mechanistic action are still elusive, and can be explored at a molecular level. Here, we employed a varied computational arsenal, including conventional and accelerated methods, to identify quinoline compounds' most likely binding sites. Our study revealed A392 and I261 as the mutations that can render RdRp resistant against quinoline compounds. In particular, for ligand 2h, mutation of A392E is the most probable mutation. The loop L1 and linker of the fingertip is recognized as a pivotal structural determinant for the stability and escape of quinoline compounds. Overall, this work demonstrates that the quinoline inhibitors bind at the template entrance channel, which is governed by conformational dynamics of interactions with loops and linker residues, and reveals structural and mechanistic insights into inhibition phenomena, for the discovery of improved antivirals.

Using Different Surface Energy Models to Assess the Interactions between Antiviral Coating Films and phi6 Model Virus.

High molecular weight chitosan (HMWCh), quaternised cellulose nanofibrils (qCNF), and their mixture showed antiviral potential in liquid phase, while this effect decreased when applied to facial masks, as studied in our recent work. To gain more insight into material antiviral activity, spin-coated thin films were prepared from each suspension (HMWCh, qCNF) and their mixture with a 1:1 ratio. To understand their mechanism of action, the interactions between these model films with various polar and nonpolar liquids and bacteriophage phi6 (in liquid phase) as a viral surrogate were studied. Surface free energy (SFE) estimates were used as a tool to evaluate the potential adhesion of different polar liquid phases to these films by contact angle measurements (CA) using the sessile drop method. The Fowkes, Owens-Wendt-Rabel-Kealble (OWRK), Wu, and van Oss-Chaudhury-Good (vOGC) mathematical models were used to estimate surface free energy and its polar and dispersive contributions, as well as the Lewis acid and Lewis base contributions. In addition, the surface tension SFT of liquids was also determined. The adhesion and cohesion forces in wetting processes were also observed. The estimated SFE of spin-coated films varied between mathematical models (26-31 mJ/m2) depending on the polarity of the solvents tested, but the correlation between models clearly indicated a significant dominance of the dispersion components that hinder wettability. The poor wettability was also supported by the fact that the cohesive forces in the liquid phase were stronger than the adhesion to the contact surface. In addition, the dispersive (hydrophobic) component dominated in the phi6 dispersion, and since this was also the case in the spin-coated films, it can be assumed that weak physical van der Waals forces (dispersion forces) and hydrophobic interactions occurred between phi6 and the polysaccharide films, resulting in the virus not being in sufficient contact with the tested material during antiviral testing of the material to be inactivated by the active coatings of the polysaccharides used. Regarding the contact killing mechanism, this is a disadvantage that can be overcome by changing the previous material surface (activation). In this way, HMWCh, qCNF, and their mixture can attach to the material surface with better adhesion, thickness, and different shape and orientation, resulting in a more dominant polar fraction of SFE and thus enabling the interactions within the polar part of phi6 dispersion.

In Silico Optimization of SARS-CoV-2 Spike Specific Nanobodies.

BACKGROUND: The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has spread worldwide, caused a global pandemic, and killed millions of people. The spike protein embedded in the viral membrane is essential for recognizing human receptors and invading host cells. Many nanobodies have been designed to block the interaction between spike and other proteins. However, the constantly emerging viral variants limit the effectiveness of these therapeutic nanobodies. Therefore, it is necessary to find a prospective antibody designing and optimization approach to deal with existing or future viral variants. METHODS: We attempted to optimize nanobody sequences based on the understanding of molecular details by using computational approaches. First, we employed a coarse-grained (CG) model to learn the energetic mechanism of the spike protein activation. Next, we analyzed the binding modes of several representative nanobodies with the spike protein and identified the key residues on their interfaces. Then, we performed saturated mutagenesis of these key residue sites and employed the CG model to calculate the binding energies. RESULTS: Based on analysis of the folding energy of the angiotensin-converting enzyme 2 (ACE2) -spike complex, we constructed a detailed free energy profile of the activation process of the spike protein which provided a clear mechanistic explanation. In addition, by analyzing the results of binding free energy changes following mutations, we determined how the mutations can improve the complementarity with the nanobodies on spike protein. Then we chose 7KSG nanobody as a template for further optimization and designed four potent nanobodies. Finally, based on the results of the single-site saturated mutagenesis in complementarity determining regions (CDRs), combinations of mutations were performed. We designed four novel, potent nanobodies, all exhibiting higher binding affinity to the spike protein than the original ones. CONCLUSIONS: These results provide a molecular basis for the interactions between spike protein and antibodies and promote the development of new specific neutralizing nanobodies.

Interfacial water molecules contribute to antibody binding to the receptor-binding domain of SARS-CoV-2 spike protein.

Antibodies that recognize the spike protein of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), especially the neutralizing antibodies, carry great hope in the treatment and final elimination of COVID-19. Driven by a synchronized global effort, thousands of antibodies against the spike protein have been identified during the past two years, with the structural information available at atomistic detail for hundreds of these antibodies. We developed an improved molecular mechanics/Poisson-Boltzmann surface area (MM/PBSA) method including explicitly treated interfacial water to calculate the binding free energy between representative antibodies and the receptor binding domain (RBD) domain of SARS-COV-2 spike proteins. We discovered that explicit treatment of water molecules located at the interface between RBD and antibody effectively improves the results for the WT and variants of concern (VOC) systems. Interfacial water molecules, together with surface and internal water molecules, behave drastically from bulk water and exert peculiar impacts on protein dynamics and energy, and thus warrant explicit treatment to complement implicit solvent models. Our results illustrate the importance of including interfacial water molecules to approach efficient and reliable prediction of binding free energy.Communicated by Ramaswamy H. Sarma.

Cu(II)(PhOMe-Salophen) Complex: Greener Pasture Biological Study, XRD/HAS Interactions, and MEP

Abstract: PhOMe-salophen (1b) (salophen is N,N-bis(salycilidene)-1,2-phenylenediamine with two tert-butyl on each ring) and Cu(II) complex with PhOMe-salophen (1c) have been synthesized and characterized using various tools, including X-ray diffraction for the Cu(II)-complex (1c, C43H52CuN2O3)). The copper complex has been obtained by Cu2+ templated approach using 1b. PhOMe-salophen (1b) has been obtained in reasonably high yield using a mixture of the Schiff-base, 1a, Pd(OAc)2, PPh3, Na2CO3, 4-methoxyphenylboronic acid in benzene. We focus in this research work on the electronic and structural properties of the Cu–Schiff base complex. The tetra-coordinate τ4 index was calculated, indicating almost a perfect square planner in agreement with X-ray diffraction results. MEP reveals the maximum positive regions in 1/-associated with the azomethine and methoxyphenyl C–H bonds with an average value of 0.03 a.u. Hirshfeld surface analysis (HSA) was also studied to highlight the significant inter-atomic contacts and their percentage contribution through 2D Fingerprint plot. In a fair comparative molecular docking study, 1b and 1c were docked together with N-[{(5-methylisoxazol-3-yl)-carbonyl}alanyl}-l-valyl]-N1-((1R,2Z)-4-(benzyloxy)-4-oxo-1-[{(3R)-2-oxopyrrolidin-3-yl}methyl]but-2-enyl)-l-leucinamide, N3 against main protease Mpro, (PDB code 7BQY) using the same parameters and conditions. Interesting here to use the free energy, in silico, molecular docking approach, which aims to rank our molecules with respect to the well-known inhibitor, N3. The binding scores of 1b, 1c, N3 are –7.8, –9.0, and –8.4 kcal/mol, respectively. These preliminary results propose that ligands deserve additional study in the context of possible remedial agents for COVID-19. © 2022, Pleiades Publishing, Ltd.

Chemoinformatics and bioinformatics by discrete mathematics and numbers: an adventure from small data to the realm of emerging big data

Currently, we are witnessing the emergence of big data in various fields including the biomedical and natural sciences. The size of chemoinformatics and bioinformatics databases is increasing every day. This gives us both challenges and opportunities. This chapter discusses the mathematical methods used in these fields both for the generation and analysis of such data. It is emphasized that proper use of robust statistical and machine learning methods in the analysis of the available big data may facilitate both hypothesis-driven and discovery-oriented research. © 2023 Elsevier Inc. All rights reserved.

Vitamin D3 and its hydroxyderivatives as promising drugs against COVID-19: a computational study.

The epidemiologic correlation between the poor prognosis of SARS-CoV-2 infection and vitamin D deficiency has been observed worldwide, however, their molecular mechanisms are not fully understood. In this study, we used combined molecular docking, molecular dynamics simulations and binding free energy analyses to investigate the potentials of vitamin D3 and its hydroxyderivatives as TMPRSS2 inhibitor and to inhibit the SARS-CoV-2 receptor binding domain (RBD) binding to angiotensin-converting enzyme 2 (ACE2), as well as to unveil molecular and structural basis of 1,25(OH)2D3 capability to inhibit ACE2 and SARS-CoV-2 RBD interactions. The results show that vitamin D3 and its hydroxyderivatives are favorable to bind active site of TMPRSS2 and the binding site(s) between ACE2 and SARS-CoV2-RBD, which indicate that vitamin D3 and its biologically active hydroxyderivatives can serve as TMPRSS2 inhibitor and can inhibit ACE2 binding of SARS-CoV-2 RBD to prevent SARS-CoV-2 entry. Interaction of 1,25(OH)2D3 with SARS-CoV-2 RBD and ACE2 resulted in the conformation and dynamical motion changes of the binding surfaces between SARS-CoV-2 RBD and ACE2 to interrupt the binding of SARS-CoV-2 RBD with ACE2. The interaction of 1,25(OH)2D3 with TMPRSS2 also caused the conformational and dynamical motion changes of TMPRSS2, which could affect TMPRSS2 to prime SARS-CoV-2 spike proteins. Our results propose that vitamin D3 and its biologically active hydroxyderivatives are promising drugs or adjuvants in the treatment of COVID-19. Communicated by Ramaswamy H. Sarma.

Design, synthesis, and molecular dynamics simulation studies of quinoline derivatives as protease inhibitors against SARS-CoV-2.

A new series of quinoline derivatives has been designed and synthesized as probable protease inhibitors (PIs) against severe acute respiratory syndrome coronavirus 2. In silico studies using DS v20.1.0.19295 software have shown that these compounds behaved as PIs while interacting at the allosteric site of target Mpro enzyme (6LU7). The designed compounds have shown promising docking results, which revealed that all compounds formed hydrogen bonds with His41, His164, Glu166, Tyr54, Asp187, and showed π-interaction with His41, the highly conserved amino acids in the target protein. Toxicity Prediction by Komputer Assisted Technology results confirmed that the compounds were found to be less toxic than the reference drug. Further, molecular dynamics simulations were performed on compound 5 and remdesivir with protease enzyme. Analysis of conformational stability, residue flexibility, compactness, hydrogen bonding, solvent accessible surface area (SASA), and binding free energy revealed comparable stability of protease:5 complex to the protease: remdesivir complex. The result of hydrogen bonding showed a large number of intermolecular hydrogen bonds formed between protein residues (Glu166 and Gln189) and ligand 5, indicating strong interaction, which validated the docking result. Further, compactness analysis, SASA and interactions like hydrogen-bonding demonstrated inhibitory properties of compound 5 similar to the existing reference drug. Thus, the designed compound 5 might act as a potential inhibitor against the protease enzyme.Communicated by Ramaswamy H. SarmaHighlightsQuinoline derivatives have been designed as protease inhibitors against SARS-CoV-2.The compounds were docked at the allosteric site of SARS-CoV-2-Mpro enzyme (PDB ID: 6LU7) to study the stability of protein-ligand complex.Docking studies indicated the stable ligand-protein complexes for all designed compounds.The Toxicity Prediction by Komputer Assisted Technology protocol in DS v20.1.0.19295 software was used to evaluate the toxicity of the designed quinoline derivatives.Molecular dynamics studies indicated the formation of stable ligand-Mpro complexes.

Computational identification of potential inhibitory compounds in Indian medicinal and aromatic plant species against major pathogenicity determinants of SARS-CoV-2.

SARS-CoV-2 (COVID-19) viral pandemic has been reported across 223 countries and territories. Globalized vaccination programs alongside administration of repurposed drugs will assumingly confer a stronger and longer individual specific immune protection. However, considering possible recurrence of the disease via new variants, a conveniently deliverable phytopharmaceutical drug might be the best option for COVID-19 treatment. In the current study, the efforts have been made to identify potential leads for inhalation therapy as nasal swabs have been reported to transfer viral load prominently. In that direction, 2363 Essential oil (EOs) compounds from Indian medicinal and aromatic plants were screened through docking analysis and potential candidates were shortlisted that can interfere with viral pathogenicity. The main protease (Mpro) of SARS-CoV-2 interacted closely with jatamansin (JM), 6,7-dehydroferruginol (FG) and beta-sitosterol (BS), while Papain-like Protease (PLpro) with friedelane-3-one (F3O) and lantadene D (LD) independently. Reduced Lantadene A (LAR) exhibited preferable interaction with RNA-dependent-RNA-polymerase (RdRp) whereas Lantadene A (LA) with RdRp and spike-glycoprotein (SG-pro) both target proteins. When compared against highest binding affinity conformations of well-known inhibitors of targets, these prioritized compounds conferred superior or comparable SARS-CoV-2 protein inhibition. Additionally, promising results were noted from pharmacokinetics prediction for all shortlisted compounds. Besides, molecular dynamics simulation for 100 ns in two replicates and binding free energy analysis revealed the stability of complexes with optimum compactness. To the best of our knowledge, the current investigation is a unique initial attempt whereby EO compounds have been computationally screened, irrespective of their known medicinal properties to fight COVID-19 infection.Communicated by Ramaswamy H. Sarma.